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Abstract 



Recently CLEO has studied the radiative decay of T into rj' and an upper limit 
for the decay has been determined. Confronting with this upper limit, most of 
theoretical predictions for the decay fails. After briefly reviewing these predic- 
tions we re-examine the decay by separating nonperturbative effect related to the 
quarkonium and that related to rj' or rj, in which the later is parameterized by dis- 
tribution amplitudes of gluons in rj' . With this factorization approach we obtain 
theoretical predictions which are in agreement with experiment. Uncertainties in 
our predictions are discussed. The possibly largest uncertainties are from relativis- 
tic corrections for J/^f and the value of the charm quark mass. We argue that the 
effect of these uncertainties can be reduced by using quarkonium masses instead 
of using quark masses. An example of the reduction is shown with an attempt to 
explain the violation of the famous 14% rule in radiative decays of charmonia. 

PACS numbers: 13.25Gv,14.40.Gx, 14.40.Aq, 12.38. Bx 
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The gluon content of 77 and 77' has been studied extensively in the literature. For example, 
recent works on the subject can be found in Radiative decays of 1 quarkonium into 77(77') 
provide an ideal place to study this subject, because the decays are mediated by gluons and there 
is no complication of interactions between light hadrons. Recently, CLEO has studied the decay 
T — ► 7 + rj' and an upper limit is determined || : 

Br(T — > 7 + V) < 1-6 x 10~ 5 (1) 

at 90% C.L.. With this result most of theoretical predictions deliver a branching ratio which is 
too large. 

The radiative decay has been studied in different approaches. In [|j] both the quarkonium and 
77(77') are taken to be nonrelativistic two-body systems, wave-functions for these bound systems 
are introduced. The obtained branching ratio in this approach with a recent compilation of a s 
is 5 — 10 x 10~ 5 [HH and is significantly larger than the upper limit. In || possible mixing 
between 77(77') and 77^ is assumed to be responsible for the decay, the branching ratio is obtained 
as 6 x 10~ 5 , which is also larger than the upper limit. In this approach it is possible to obtain 
Br(T — > 7 + 77') pa (1 ~ 3) x 10~ 5 close the upper bound 0. 

The corresponding decay of J/ty has been studied by saturating a suitable sum rule with J 
resonance and it has been shown that the decay is controlled by the C^a(1) anormaly J/J. The 
result of this study can be re-written in the form: 

r(j/» -> 7 + n') -^.{i-rgLf. i|(o|^C^hrt|T-' W * - «+«-), (2) 

where G a,tiu is the field strength tensor of gluon and G a ^ v = jE^ ual3 G a a/3 . In the above result we 
have neglected the binding energy of J/^f and taken Mj/^ = 2m c , m c is the pole mass of c-quark. 
If one can generalize the approach for the T decays, one can obtain the ratio: 
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Using the experimental results for the widths in the bracket we obtain 
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The branching ratio of J/\I/ has been measured and its value is (4.31 ± 0.3) x 10~ 3 . We take the 
quark masses as m& = My/2 pa 5GeV and m c = Mj/^/2 ~ 1.5GeV and obtain Br(T ^7 + 77') pa 
3.1 x 10 -4 , which is too large for the upper limit. However, the generalization of Eq.(2) to T may 
not be correct. In the spirit of the approach the emitted gluons, which are converted into 77', are 
soft, while in T decay the gluons are definitely hard. Employing multipole expansion for the soft 
gluons one is also able to predict the decay of J/ty |§ 
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From above discussions one may conclude that the predictions based on QCD-inspired models 
or on sum rule are not compatible with the upper limit, or not consistent. It should be also noted 
that phenomenological models can have compatible predictions. In an extended vector-dominance 
model one indeed finds the branching ratio from 5.3 x 10~ 7 to 2.5 x 10~ 6 ||, but this model has 
no direct relation to QCD as the fundamental theory of strong interaction. 

Decays of quarkonia were intensively studied in eighties. Now our understanding of QCD has 
been greatly improved, a restudy of these decays is necessary to explain new experimental results 
like the upper bound in Eq.(l). On the other hand, there is a large data sample with 5 x 10 7 
J/ip events collected with BES HlQ| , a date sample with several billions J/\I/ events is planed to be 
collected with the proposed BES III at BEPC II and with CLEO-C at modified Cornell Electron 
Storage Ring (CESR) Furthermore, about 4 fb" 1 bb resonance data are planned to be 



taken at CLEO III in the year prior to conversion to low energy operation (CLEO-C) JO]. These 
data samples of quarkonia will allow us to study the decays, which have been observed before, 
with more accuracy, and also those decays which have been not observed. Therefore, experimental 
activities will bring more information about these decays and may also lead to new discoveries, 
e.g., discovery of glueball. In this brief report we present an approach based on QCD factorization 
to explain the experimental result from CLEO 0. This approach was used for the radiative decay 
into the tensor meson fa [12| . 



We consider the heavy quark limit, i.e., m c — > oo, m& — > oo. In the limit, a quarkonium 
system, taking T as an example, can be taken as a bound state of b- and 6-quark which move 
with a small velocity v, hence an expansion in v can be employed, nonrelativistic QCD(NRQCD) 
can be used to describe the nonperturbative effect related to T The decay can be thought 



as the following: the quarkonium will be annihilated into a real photon and gluons, the gluons 
will be subsequently converted into the meson rf . Also in the limit, the meson rj' has a large 
momentum, this enables an expansion in twist to characterize the gluonic conversion into rf, the 
conversion is then described by a set of distribution amplitudes of gluons. The large momentum 
of rf requires that the gluons should be hard, hence the emission of the gluons can be handled by 
perturbative theory. The above discussion implies that we may factorize the decay amplitude into 
three parts: the first part consists of matrix elements of NRQCD representing the nonperturbative 
effect related to T, the second part consists of some distribution amplitudes, which are for the 
gluonic conversion into rf , the third part consists of some coefficients, which can be calculated 
with perturbative theory for the 66-pair annihilated into gluons and a real photon. In this report 
we show that the contribution of twist- 2 operators are suppressed by rnii. This indicates that a 
complete QCD-analysis should include contributions from twist-4 operators. However, without 
such a complete analysis one still can make some predictions like the branching ratio given in 
Eq.(l). 

We consider the decay of T 

T->7(«) +»/(*), (5) 

where the momenta are given in the brackets. We take a light-cone coordinate system, in which 
the momentum k of rf is k^ = (k + , k~ , 0, 0). We consider the contribution from emission of two 
gluons, and assume a factorization can be performed. Then the S-matrix can be written as: 



2 



(1V'\S\T) = -i\eQ h g% J d'xd'yd'zd'x, d 4 yi e^ z ( V '\G^x)G b u (y)\0) 

(Olb^MyJlT) ■ M?r' ab (x,y,x u y u z), (6) 

where Mj^ p,ab (x, y, x±, yi, z) is a known function, i and j stand for Dirac- and color indices, a 
and 6 is the color of gluon field, b(x) stands for the Dirac field of 6-quark, e* is the polarization 
vector of the photon, Q b is the charge fraction of the 6-quark in unit e. The above equation can 
be generalized to emission of arbitrary number of gluons. Using the fact that 6-or 6-quark moves 
with a small velocity v, the matrix element with the Dirac fields can be expanded in v. We obtain: 

(0|^)My)|T> = ~(P + ^P.)Mx^\T)e^^ + 0(v 2 ), (7) 

where i s the NRQCD field for 6(6) quark and 

P ± = (l± 7 °)/2, 

p" = (m 6 , 0,0,0), (8) 

where m b is the pole-mass of the 6-quark. In Eq.(7) we do not count the power of v for quark 
fields because this power is same for every term in the expansion of Eq.(7). With this in mind the 
leading order of the matrix element is then O(v ), we will neglect the contribution from higher 
orders and the momentum of T is then approximated by 2p. It should be noted that effects at 
higher order of v can be added with the expansion in Eq.(7). 

For the matrix element with gluon fields we observe that the x-dependence of the matrix 
element is controlled by different scales: the ^'-dependence is controlled by k + , while the x + - 
and a^-dependence is controlled by the scale Aqcd or k~, which are small in comparison with 
k + . Because of these small scales we can expand the matrix element in x + and in xt- With this 
expansion we obtain the result for Fourier transformed matrix element: 

j dx'e-^WlG^G^mO) = l -5 a \2n)%q^)5\q lT ) ■ ^-1— ^ F,,^) 

+ •••, 1 1 (9) 

with 

= e^ a %n p , l» = (1,0,0,0), n^ = (0,1,0,0), qf = x,k + , 
F ri ,(x 1 ) = ^ J dx-e- lXlk+x ~( V \k)\G a > + ^x-)G a > + »(0)\0) eilu , (10) 

where e^ al3 is totally anti-symmetric with e 0123 = 1. F Tj >(x 1 ) is the distribution amplitude char- 
acterizing the conversion of two gluons into rj' and it is defined with twist-2 operators in the 
light-cone gauge. In other gauges a gauge link should be supplied in Eq.(10) to maintain the 
gauge invariance. It should be noted that there is no simple relation between F v /(x) and the 
gluonic matrix element in Eq. (2). The • • • in Eq.(9) stands for contributions from higher twist. 
The next-to-leading twist is 4. With the above results we obtain the ^-matrix element with the 
twist-2 contribution in the limit m v > — > 0: 
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(V/W) = ^eQ^ s 2 (2vr) 4 ( 5 4 (2p - k - g)e;<0| X Vty|T) 

Xi(l — X\ 




<> \ dx K JTZr^Y2 F n'M ■ ( 1 + 0{^) ) . (11) 



It shows that the twist-2 contribution is suppressed by m 2 ,. In the twist-expansion the light 
hadron mass m v > should be taken as a small scale as Aqcd, hence the contribution is proportional 
to A.qcd- This implies that a complete analysis at the leading order should include not only this 
contribution but also twist-4 contributions, in which one needs to consider the contributions from 
emission of 2, 3 and 4 gluons. This is too complicated to be done here. However, without a 
complete analysis we can always write the result of a complete analysis as 

(j V '\S\T) = ^eQ b gl{2^5\2p - k - g) e ;<0| X Vty|T) e^-L . ^ (12) 

where the parameter g v / has a dimension 3 in mass. This parameter is a sum of the twist-2 
contribution in the second line of Eq.(ll) and the twist-4 contributions which need to be analyzed. 
The parameter characterized the conversion of gluons into rf and it does not depend on properties 
of T. The origin of the factor m 4 = ml ■ ml is: one ml comes from the perturbative part, 
another ml reflects the fact that the contribution of twist-2 operators is proportional to m 2 , and 
contributions of twist-4 are proportional to Aq CD . It is interesting to note that this power behavior 
is also obtained in JTJJ, in contrast, it is also pointed out in |l4j that this behavior holds by taking 



m c and m& as light quark masses, in the heavy quark limit this behavior does not hold. 

The above result may also be generalized for the decay of J/^ . One may have a doubt that 
the twist expansion may not be applicable for the J/^> decay, because m c is not large enough. 
The twist expansion means a collinear expansion of momenta of partons in rj', components of 
these momenta have the order of (0(k + ), 0((k~), O(Aqcd), 0(Aqcd))- Hence the expansion 
parameters are 

k~ ml< Aqcd , s 

F=5fr 01 ' -f^ - 2 - (i3) 

where we have taken Aqcd ~ 400MeV. These estimations show that the twist expansion is also 
a good approximation for the J/^f decay. There are also other possible large uncertainties, due 
to effects from higher orders of v and a s . These can be eliminated partly by building the ratio 
r(J/ 1 J r (T) — > 77/)/r(J/ 1 i r (T) — > light hadrons). Theoretical prediction for this ratio will have less 
uncertainties than the width, because corrections from higher orders of v and a s are cancelled 
at certain level. Using experimental data for Br(J/ 1 i r (T) — > light hadrons) we can predict the 
branching ratio. With this consideration we rewrite the ratio defined in Eq.(3) as: 

Br(T — >• light hadrons) 
~ Br(J/$ -»• light hadrons) ' ^'^ 

r(T — > jt]') /r(T — > light hadrons) Ql m c ft s(m c ) 
Tr) ' ~ r( J/$ -»• 7r/)/r(J/# ~> light hadrons) ^ Q 2 c m 6 b ' a s (m ft )' ^ ' 
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where leading order results for the decay widths are used for r v i. Using the experimental results 
for the branching ratios of decays into light hadrons, we obtain: 

Br(T -*_2 + rf) Q%m 6 c a s (m c ) 

Br( - 7 + rf) ~ i dl ' Q2 m 6 • as(mfc) • W 

This is the result at the leading order of A, where A is Aqqd or m^, and the dependence of the 
renormalization scale in gluonic distribution amplitudes is neglected. The dependence may be 
extracted from the study in [[DJ. By taking a s (m c ) ~ 0.3 and a s (rrib) ~ 0.18 we obtain: 

R v , « 3.9 x 10" 4 . (16) 

With the experimental value of Br( J/\]/ — > 7 + 7/) we obtain the branching ratio: 

Br(T ^7 + 7/) « 1.7 x 10~ 6 . (17) 

This value is much smaller than the values obtained with other approaches and it is in consistency 
with the upper limit. Similarly we also obtain 

Br(T -> 7 + 77) « 3.3 x 10" 7 . (18) 

It should be emphasized that our results obtained in the above equations are not based on any 
model, corrections to these results can be systematically added in the framework of QCD. The 
possibly largest uncertainties in our results are from relativistic corrections for <//*& decays and 
from the uncertainty of the value of the charm quark mass, each of them can be at the level of 50%. 
Taking these into account, our prediction in Eq.(17) and (18) can be close to the experimental 
bound in Eq.(l). However, these largest uncertainties may be reduced by using hadron masses, 



i.e., using 2m c = Mj/\&. This possibility is based on the result for relativistic correction in |L6 
and on the observation that the violation of the famous 14% rule may be reduced in this way. 
If one analyzes the correction at the next-to leading order of v for decays of 1 quarkonia, one 
obtains that the correction is proportional to a NRQCD matrix element defined in ||13|| . This 
matrix element represents a relativistic correction. In ]Tj| it is shown that this matrix element 
is proportional to the binding energy, i.e., to Mj/^ — 2m c for J/ty and to Afy — 2m c for ip' 
respectively. If we use 2m c = Mj/^, for J/ty decays and 2m c = M^> for ip' decays respectively, 
the relativistic correction disappears formally, but it is actually included by using hadron masses. 
However, it should be noted that this should be regarded as a phenomenological estimation, a 
detailed analysis and an precise determination of quark masses is needed to study the correction 
in a consistent way. 

The famous 14% rule is derived simply by taking leading order results for decays. In our case 
we have: 

Br(-?// — > 777') Br(ib' e + e~) , , 

111 - ' -- 0.147 ±0.023, (19) 



Br(J/^ -> 7 r/') Br( J/tf -> e+ e - 

where the number is estimated with experimental results of leptonic decay widths. This result is 
theoretically expected not only for radiative decays into any light hadron, but also for hadronic 
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decays, this is the so-called 14% rule. However, this rule is significantly violated, one of the 
violations is the well known pn puzzle. A possible explanation and useful references can be found 
in ]T7|]. The experimental result made by BES indicates that the rule is also violated in our 

Br(^// — > 777') 



case: 



Br( J/m -> 77/) 



0.036 ±0.009. (20) 



This value is only fourth of the expected. It should be noted that corrections from higher orders 
of a s is cancelled in the ratios in Eq.(18), the theoretical uncertainties come from effects of higher 
orders in v in Eq.(7) and higher twists. In the case of 1 quarkonia, the correction from the 
next-to-leading order of v is the relativistic correction, whose effect is expected to be significant 
for charmonia. As discussed before, this correction may be estimated by replacing m c with the 
half of the mass of quarkonium, i.e., we use 2m c = Mj/^ for the J/ty decays and 2m c = for 
if)' decays. With this replacement and with our result in Eq.(12), the ratio in Eq.(19) is modified 
as: 



Br<y^7?/) Mj/* BrO' ^ e + e 



Br( J/V -»■ 77/') M|, Br( J/* -> e+e 



0.0512 ±0.0080, (21) 



This result shows that the relativistic correction is indeed significant. With the replacement 
the predicted ratio is much closer to the experimental result than that of the 14% rule and the 
two largest uncertainties are reduced in the prediction. However this is a naive estimation, a 
detailed study is needed and is in progress |19| . With this case we can expect that the two 
largest uncertainties are also reduced in our predictions in Eq.(17) and (18) because we have used 
2mfe = M-f and 2m c = Mj/^. 

It is also interesting to look at decays into pir. In this decay one of the final hadrons is 
produced at the level of twist-2, another is at the level of twist-3 [pOfl . With this fact and with 
the replacement the rule is modified as: 



Br^'^pvr) Mj^ Bi(ip' 



e T e 



Qp* = ttTtTt — — ~\ = * A T~i / t / t — — ~t = 0.036 ± 0.006. (22) 

^ p Bi(J/m pn) M$, Br(J/^ -> e+e") V ; 

With the modification the rule is changed significantly. The above result also holds for decays 
into K*K. Although the ratio is reduced, but it is still in conflict with experimental results. In 
21| it is found that Q pn < 0.006 and Q k *+k- < 0.64. Recent data from BES gives Q pn < 0.0023 



and Qk*°k° — 0.017 ± 0.006 p2| . However, the predictions are closer to experimental than that 
in Eq.(19). One should also keep in mind that these decays are more complicated than radiative 
decays discussed before, because the final state consists of two light hadrons. 

To summarize: We have presented a QCD-factorization approach for radiative decays of 1 
quarkonium into i](r]'), the result is in consistency with the experimental result made by CLEO. 
On the other hand, most of theoretical results is not compatible with the upper limit. A possible 
explanation for the violation of the 14% in our case is given. With this explanation we show that 
effect of relativistic corrections and that due to uncertainty of the quark mass can be reduced by 
using quarkonium masses and uncertainties in our predictions may be not so large as those usually 
expected. 
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